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Abstract. This article is devoted to results on 7C~-7T~-Bose-Einstein correlations in central Pb+Pb 
collisions measured by the NA49 experiment at the CERN SPS. Rapidity as well as transverse 
momentum dependences of the correlation lengths will be shown for collisions at 20A, 30A, 40A, 
80A, and 158A GeV beam energy. Only a weak energy dependence of the radii is observed at SPS 
energies. The ^-dependence of the correlation lengths as well as the single particle m,-spectra will 
be compared to model calculations. The rapidity dependence is analysed in a range of 2.5 units of 
rapidity starting at the center of mass rapidity at each beam energy. The correlation lengths measured 
in the longitudinally comoving system show only a weak dependence on rapidity. 
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1. INTRODUCTION 

The measurement of correlations of identical bosons in heavy ion collisions provides 
a unique tool to investigate the space time evolution of the particle emitting source. 
Bose-Einstein correlations are observed as an enhancement of the yield of pairs of 
particles with small relative momenta. Measurements of the range and strength of the 
correlations in momentum space allow to derive the extension of the source in coordinate 
space. Due to space momentum correlations in expanding sources the correlations do 
not reflect the whole extensions of the source. In such a scenario, the study of the 
correlations in different regions of phase space helps to understand the evolution of 
the source. While the dependence of the correlation lengths on the mean transverse 
momentum k t = ^\pt,i + Pt.i\ of the pairs reflects the transverse expansion dynamics 

of the source the dependence on the pair rapidity Y = j log ( ^^^z^—z^- j > which is 
measured in the center of mass system, should shed light on the profile of the source in 
longitudinal direction. The large acceptance of the NA49 experiment allows us to obtain 
a comprehensive picture of the dynamical evolution of the source. 

The article is organized as follows: A brief survey of the experiment, the construction 
of the correlation function and the fit method are presented in section 2. In section 3, 
crucial systematic uncertainties in the analysis due to detector effects are discussed. The 
results on the k t - as well as the Y- dependence of the correlation lengths are presented in 
section 4 and compared to model calculations. 

2. EXPERIMENTAL SETUP AND ANALYSIS 

NA49 [1] is a fixed target experiment located at the CERN SPS comprising four large- 
volume Time Projection Chambers (TPC), two of which are located inside the magnetic 
field of two superconducting dipole magnets (Figure 1). The TPCs are read out at 90 
(MTPC) and 72 (VTPC) pad rows resulting in a very good determination of the momen- 
tum of the traversing particles. A zero degree calorimeter at the downstream end of the 
experiment is used to trigger on the centrality of the collisions. The data presented here 
correspond to the 7.2% most central events for data samples taken at 20A, 30A, 40A, 



13 m 




TARGET 



MTPC-R /7 1 
TIME PROJECTION CHAMBERS 



CALORIMETER 



FIGURE 1. NA49 detector setup. 



80A, and 158A GeV beam energy. By scaling the magnetic field it was possible to obtain 
a similar coverage of phase space relative to the center of mass rapidity for the different 
beam energies. Measuring the specific energy loss dE/dx of charged particles in the gas 
of the TPCs with a resolution of 3-4% allows particle identification. However, due to 
ambiguities in particle identification by specific energy loss measurements in certain re- 
gions of phase space, negative hadrons rather than identified negative pions are studied 
in this analysis. 

The correlation function is constructed as the ratio of a distribution of the momentum 
difference of pairs from the same event (signal) and a mixed event background distribu- 
tion (background). Following the approach of Pratt and Bertsch [2, 4] the momentum 
difference is decomposed into a component parallel to the beam axis qi ong and two 
components in the transverse plane q out and q S id e with q out defined parallel, and q si( i e 
perpendicular to k t . The correlation function is parameterised by a Gaussian function 

2 2 2 2 2 2 2 

C2{q)BP = 1 + A ■ exp(-R out q < out -R side q side -Rl on g1long ~ ^Routlongloutqiong) (1) 

and the parameters R out , R s ide, Riong, Routiong, and A are determined by a fit to the 
measured correlation function. The Coulomb repulsion of the particles is accounted 
for by weighting the theoretical Bose-Einstein correlation function C2(q)BP by a fac- 
tor F(qi m „<r>) [3] in the fit procedure. The weight is determined by the invariant mo- 
mentum difference qi m of the pair and the mean pair separation <r> of the particles in 
the source. According to [3] this quantity can be derived from the extracted radii. We 
therefore determine the source parameters as well as the value of <r> in an iterative 
fit-procedure. For following fit function was used: 



C 2 (q) f = n{p- (C 2 (q) B p • F(q, < r >)) + (1 - p)} . 



(2) 



The contamination of the sample with pairs of non-identical particles and pairs of pions 
from long lived resonances or weak decays is accounted for by a purity factor p which 
is determined by a VENUS/GEANT simulation. The fit parameter n is introduced to 
account for the different statistics in signal and background distributions. Beside the 
uncertainties which arise due to the construction of the correlation function and the fit 
formalism there are further detector related effects which will be discussed in the next 
section. 



3. SYSTEMATIC STUDIES 



3.1. Two track resolution 

The momentum difference of a pair is closely related to the distance of the tracks 
traversing the detector. Bose-Einstein correlations are restricted to a narrow window in 
momentum difference, hence it is crucial to understand the two track resolution of the 
detector. The overlap of charge clusters induced at the pad planes of the TPCs can lead 
to an assignment of points to the wrong track or in an extreme case to complete merging 
of two tracks. In this case, a pair with small relative momentum will be lost in the signal 
distribution and the observable Bose-Einstein enhancement will be reduced. To study 
the impact of the limited two track resolution on the extracted radii, the distance of the 
tracks was measured at each pad row where both tracks lie in the sensitive volume of the 
TPCs. Starting from the downstream end of the TPCs the distance of closest approach 
idea) of two tracks after a given number of passed pad rows n rows was determined. Pairs 
with a dca smaller than a given cut value dca cllt were rejected both from the signal and 
the background distribution. The impact of a variation of the two parameters dca cut and 
n rows on the correlation function is shown in Figure 2. Requiring only small values of 
dca cut and n rows , tracks can approach each other very closely over a considerable part of 
the track length in the TPCs. In this case, track merging effects can lead to a significant 
loss of pairs in the signal. This effect is very pronounced at high transverse momenta 
and shows up in an undershoot of the projection of the correlation function onto q out 
(Figure 2a). Increasing the cut parameters, the influence of merging effects is reduced, 
the undershoot of the correlation function vanishes and the extracted radii vary only by 
less than 0.2 fm (Figure 2b). This can serve as an estimate of the systematic error on the 
radii due to the specific treatment of the two track inefficiencies in this analysis. 

Further systematic errors on the radii arise due to uncertainties concerning the treatment 
of the Coulomb interaction (which can significantly influence the parameters R out and 
A), the missing particle identification, the finite momentum resolution of the detector and 



a) 158AGeV k t =(0.3-0.4)GeV/c Y = 0.0-0.5 b) R ut 
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FIGURE 2. a) Impact of two track resolution inefficiencies on the shape of the correlation function and 
b) fit results for R ou1 for different combinations of the cut parameters dca cu t and n rows . 



the normalisation of the correlation function. The overall systematic error on the radii 
is not specified for each k t -Y -bin but estimated to be smaller than 1 fm for all extracted 
radii. 

4. RESULTS 
4.1. The /^-dependence 

Figure 3 presents the ^-dependence of the radii at midrapidity. Radius and X param- 
eters were obtained from fits of (2) to the correlation function in bins of Y and k t . Since 
the influence of the finite momentum resolution was found to be small, no corrections 
were applied for this effect. The bin width in k t was chosen to be 0.1 GeV/c and in- 
creased to 0.2 GeV/c for the last bin ((0.4-0.6) GeV/c) to obtain sufficient statistics at 
all energies. A strong decrease of Ri ong with k t is observed. R out is slightly larger than 
R s ide for all energies indicating a finite duration of particle emission [4]. Also shown 
in this picture is a fit of the ^-dependence of the radii according to a blast wave pa- 
rameterisation [5] of the source. In this model the source is treated as boost invariant in 
the longitudinal direction. In the transverse direction a box-shaped density profile and 
a linearly increasing flow profile is assumed. Space momentum correlations induced by 
flow reduce the measured correlation lengths. This effect is partly compensated in case 
of a superimposed thermal velocity field. Therefore ambiguities arise in the two model 
parameters temperature and flow, which can not be resolved by only analysing the k t - 
dependence of the radii. To resolve these ambiguities the single particle p t - spectra of 
protons and negatively charged pions measured by NA49 were fitted to a parameterisa- 
tion derived from the same model. The lines in Figure 3 correspond to a combined fit 
to the radii and the particle p r spectra. As expected from the weak energy dependence 
of the radii only small variations of the extracted source parameters were observed. The 
extracted parameters were the temperature T, the maximum transverse flow rapidity p, 
the transverse geometrical radius R, the emission time x and the emission duration At. 
The fit results are inserted in Figure 3. The temperature T slightly increases with the 
beam energy, transverse flow and geometrical radius stay approximately constant over 
the observed energy range. For the emission time we obtain values of 5.4 to 6.8 fm/c. 
The fit slightly overpredicts R si( i e at high k t but still results in a finite emission duration 
of 2.2-3.2 fm/c. The fit might be further constrained by adding more particle spectra or 
by including contributions from resonance decays in the model. 



4.2. Y -dependence 

The model described in section 4.1 is only applicable in case of a longitudinally 
boost invariant source. Under such conditions it is expected that the cross term R ou tlong 
vanishes [6]. Considering the systematic error of 1 fm this condition is fulfilled to 
good approximation at midrapidity. In Figure 4 the rapidity dependence of Ri ong , R s ide> 
Rout, an d R mlong is shown at £ f =(0.0-0.1) GeV/c for the different beam energies. In 
[7] the impact of a non-boost invariant expansion on the parameter R ou tlong is studied. 
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FIGURE 3. The ^-dependence of R s a e , Rout, and Ri ong at midrapidity (0.0 < Y < 0.5) for the different 
data sets (dots). The lines correspond to a combined fit of the blast wave model to the radii and the single 
particle p r -spectra. 



An increase of R ou tlong with increasing rapidity is predicted due to the decrease of the 
inclusive pion yields and is in agreement with the results for all energies. However 
a change in the longitudinal expansion dynamics which is indicated by the change 
in Routlong is not reflected in the rapidity dependence of the other observables. R s id e , 
which is supposed to determine the geometrical size of the source [8] does not change 
significantly with rapidity. R out is approximately constant over the investigated rapidity 
region. Only slight changes in Ri ong are observed. These are even less pronounced at 
higher transverse momenta. 

In summary, a distinct energy dependence of the radii is not observed even though there 
is a dramatic change in the energy dependence of other hadronic observables like e.g. 
the kaon to pion ratio [9]. Furthermore the radii do not show a pronounced rapidity 
dependence, in contrast to the particle yields which decrease strongly with rapidity. 
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FIGURE 4. Rapidity dependence of R s id e , Rout, Rfong> an d Routlong at & f =(0.0-0. 1) GeV/c for the different 
beam energies. Shown are as well the midrapidity values Y m j ( [ for the different beam energies. 



Hungarian Scientific Research Foundation (T032648, T032293, T043514), the Hungar- 
ian National Science Foundation, OTKA, (F034707), the Polish-German Foundation, 
the Korea Research Foundation Grant (KRF-2003-070-C00015) and the Bulgarian 
National Science Fund (Fh-09/05). 



REFERENCES 

1. S. Afanasiev et al., NIM A430 (1999) 210-244 

2. S. Pratt, Phys. Rev. D 33 (1986) 1314-1327 

3. Yu. M. Sinyukov et al., Phys. Lett. B 432 (1998) 248-257 

4. G. F. Bertsch, Nucl.Phys. A 498 (1989) 173c-180c 

5. F. Retriere, M. Lisa, Phys. Rev. C 70 044907 (2004) 

6. S. Chapman, P. Scotto und U. Heinz, Phys.Rev.Lett 74 (1995) 4400-4403 

7. S. Chapman, P. Scotto und U. Heinz, Nucl.Phys. A 590 (1995) 449c-452c 

8. S. Chapman, J. Nix and U. Heinz, Phys. Rev. C 52 (1995) 2694-2703 

9. M. Gazdzicki, J. Phys. G: Nucl. Part. Phys. 30 (2004) 701-708 



